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ABSTRACT 

Some general properties of the advective- acoustic instability are described and understood using 
a toy model which is simple enough to allow for analytical estimates of the eigenfrequencies. The 
essential ingredients of this model, in the unperturbed regime, are a stationary shock and a subsonic 
region of deceleration. For the sake of analytical simplicity, the 2D unperturbed flow is parallel 
and the deceleration is produced adiabatically by an external potential. The instability mechanism is 
determined unambiguously as the consequence of a cycle between advected and acoustic perturbations. 
The purely acoustic cycle, considered alone, is proven to be stable in this flow. Its contribution to 
the instability can be either constructive or destructive. A frequency cut-off is associated to the 
advection time through the region of deceleration. This cut-off frequency explains why the instability 
favours eigenmodes with a low frequency and a large horizontal wavelength. The relation between 
the instability occurring in this highly simplified toy model and the properties of SASI observed in 
the numerical simulations of stellar core-collapse is discussed. This simple set up is proposed as a 
benchmark test to evaluate the accuracy, in the linear regime, of numerical simulations involving this 
instability. We illustrate such benchmark simulations in a companion paper. 
Subject headings: accretion - hydrodynamics - instabilities - shock waves - supernovae 



1. INTRODUCTION 

In the scenario proposed by Bethe & Wilson (1985) 
for core collapse supernovae, the success of the explo- 
sion depends on the efhciency of energy deposition by 
neutrinos below the stalled accretion shock, during the 
first second after core bounce. Numerical simulations 
revealed that this mechanism is inefficient in spherical 
symmetry (Liebendorfer et al. 2001). Hydrodynamical 
instabilities may play an important role by breaking the 
spherical symmetry and helping the revival of the stalled 
shock. Indeed, multidimensional simulations allowing for 
transverse motions, induced by neutrino-driven convec- 
tion in the gain region, approached the explosion thresh- 
old (Burrows et al. 1995, Janka & Miiller 1996), al- 
though this effect did not seem sufficient (Buras et al. 
2003). The discovery of another hydrodynamical insta- 
bility, named the Standing Accretion Shock Instability 
(SASI) by Blondin et al. (2003), opened new perspectives 
which seem very promising for the revival of the shock. 
The recent simulations of Marek & Janka (2007) showed 
a successful explosion of a ISM© progenitor where neu- 
trino energy deposition is efficient enough owing to the 
effect of SASI, which lengthens the time spent by the 
postshock gas in the gain region (Murphy & Burrows 
2008). SASI is also the starting point of the acoustic 
mechanism found by Burrows et al. (2006, 2007), where 
the excitation of g-mode oscillations inside the proto- 
neutron star triggers the emission of acoustic waves (see 
also Weinberg & Quataert 2008). Besides the question 
of the explosion mechanism, the development of SASI 
seems to have important consequences on the birth con- 
ditions of the neutron star, its kick (Scheck et al. 2004, 
2006) and also its spin (Blondin & Mezzacappa 2007, Ya- 
masaki & Foglizzo 2008). 
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In view of the spectacular possible consequences of SASI, 
a fundamental understanding of its mechanism is de- 
sired but still a source of debate. It is definitely distinct 
from convection since it can take place even without any 
source of heating (Blondin et al. 2003). The "advective- 
acoustic cycle" (Foglizzo & Tagger 2000, Foglizzo 2001, 
Foglizzo 2002) was recognized as the driving mechanism 
by several authors (Blondin et al. 2003, Burrows et al. 
2006, Ohnishi et al. 2006, Foglizzo et al. 2007, Scheck et 
al. 2008, Yamasaki & Foglizzo 2008), but some alternate 
interpretations were proposed (Blondin & Mezzacappa 
2006, Blondin & Shaw 2007). It should be noted that 
the analytical arguments of Laming (2007) supporting 
the existence of a purely acoustic mechanism contained 
some errors, and the corrected calculation favours the 
advective-acoustic mechanism (Laming 2008, in press). 
Part of the difficulty in recognizing the advective-acoustic 
mechanism in numerical simulations, even in the simpli- 
fied set up proposed by Blondin et al. (2003), comes from 
the lack of simple reference models where its properties 
would be fully understood. The present work aims at 
providing such a reference, as simple as possible. This 
reference serves two purposes. It can help us build our 
physical intuition about the advective-acoustic coupling 
responsible for an unstable advective-acoustic cycle. It 
can also be used as a benchmark test to evaluate the ac- 
curacy of multidimensional numerical simulations in the 
linear phase of the instability (Sato et al. 2008, hereafter 
paper II). 

Although some of the equations used in the present work 
have already been described in past publications dealing 
with more complex flows (e.g. Yamasaki & Foglizzo 2008 
and references therein) , we choose to explain them again 
in the simpler framework considered here. Analytical cal- 
culations are made possible by three main simplifications: 
(i) the flow is 2D planar, (ii) the region of deceleration is 
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Fig. 1. — Schematic view of the toy model. Entropy /vorticity 
perturbations (circular arrows) are advected downward with the 
flow, and coupled to acoustic ones (wavy arrows) in the inhomoge- 
neous region near z\j . The linear coupling between these pertur- 
bations is described by coupling coeflicients Qsh> Svi ^slu ^Vi 
defined in Sect. 3 



localized spatially at a distance from the shock, and (iii) 
the deceleration is produced adiabatically. The simplic- 
ity of the present set up allows for an exact calculation, 
even at low frequency. The physical insight in the mech- 
anism at work can be used as a reference, in order to 
analyze more realistic situations. 

The paper is organized as follows: in Sect. 2, the toy 
model consisting of a stationary shock in an external 
step-like potential is described, and its eigenspectrum is 
computed as an illustrative example. Sect. 3 explains the 
method used to prove unambiguously the mechanism at 
work, beyond the determination of the eigenfrequencies. 
This method requires the determination of the efficiency 
of the advective-acoustic coupling at the shock, in Sect. 4, 
and in the region of acoustic feedback in Sect. 5. These 
local wave coupling efficiencies are combined into a global 
cycle whose growth rate is computed in Sect. 6 and com- 
pared to the direct eigenspectrum of Sect. 2. The impor- 
tance of the size of the deceleration region is emphasized 
in both Sect. 5 and 6. The relevance of these results 
to the problem of core-collapse is discussed in Sect. 7. 
Results are summarized in Sect. 8. Lengthy analytical 
derivations are described in Appendices for the sake of 
clarity. 

2. DESCRIPTION OF A TOY MODEL 

2.1. Stationary flow 

The essential ingredient of the advective-acoustic in- 
stability is the interplay of advected and acoustic per- 



turbations in a subsonic cavity (e.g. the sketch drawn in 
Fig. 1 of Scheck et al. 2008). For the sake of simplic- 
ity, we build a toy model where the advective-acoustic 
cycle takes place without the complications associated 
with the spherical geometry, and without the difficulties 
relative to non adiabatic processes. The discussion of 
these difficulties is postponed to Sect. 7. We choose the 
2D parallel flow of an ideal gas in a Cartesian geometry, 
passing through a shock and adiabatically decelerated 
by an external potential. The density is uniform in the 
X direction and the gas flows along the z direction with 
a negative velocity. The initial position of the stationary 
shock is noted Zgh and the region of deceleration centered 
on zv (Fig. [T|) . The distance from the shock to the center 
of the deceleration region is noted H = Zsh — z-^ ■ 

In the adiabatic approximation, the equations describ- 
ing the stationary flow of ideal gas, with a density p and 
a velocity v in the external potential $(2;), are the con- 
servation of mass flux and the Euler equation: 
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where 7 is the adiabatic index, and the sound speed c is 
related to the pressure p by = 7p//5. The only gradi- 
ents in this stationary flow are produced by the external 
potential ^{z), localized in the vicinity of the coordinate 
zv over a region of size H-^: 
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At a distance exceeding ^ 3-ffv from z^, the flow is ap- 
proximately uniform (1 — tanh(6) ^ 1.2 x 10~^). Quanti- 
ties in the uniform subsonic regions upstream and down- 
stream of the potential jump are denoted by the sub- 
scripts "in" and "out" respectively, while the subscripts 
"1" and "sh" denote quantities immediately ahead and 
after the stationary shock respectively. The entropy S is 
defined by S" = \^og{{p/psh)/{pl PshV)] /{l - 1), and the 
Mach number is defined as positive: M. = —v/c. The 
only dimensionless parameters of the stationary flow are 
the adiabatic index 7, the incident Mach number M.i, 
the relative size of the coupling region / H and the 
adiabatic heating parameter Cq^^/c^j (or the deceleration 
parameter Wout/i'in)- Note that 7 = 4/8 and Mi — 
throughout the paper. The potential jump A$, the 
sound speed, velocity, and Mach number jumps are re- 
lated to the parameters of the flow, through conservation 
equations, as follows: 
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2.2. Linear perturbations 

The ID stationary flow is perturbed in the plane [x, z), 
with periodic boundary conditions in the direction x. 
For a given horizontal size L,j. of the computation do- 
main, the wavenumber is restricted to discrete values 
associated to the number rix of horizontal wavelengths, 
kx = 2i:nx/ Lx- The incoming supersonic flow is not per- 
turbed. For the sake of physical simplicity, the lower 
boundary condition is defined by the free leakage of per- 
turbations, i.e. the absence of an acoustic flux propa- 
gating upward from z < Z'^ — 3i?v- Instead of repeating 
the tedious derivation of the differential system governing 
the evolution of linear perturbations, we deduce it from 
previous studies such as Foglizzo et al. (2006, hereafter 
FSJ06) or Yamasaki & Foglizzo (2008), by neglecting non 
adiabatic effects or rotation. The three functions chosen 
to describe the evolution of perturbations are noted 5f, 
Sh and SS, where SS is the entropy perturbation, df is 
the the perturbation of the energy density associated to 
the Bernoulli invariant, and Sh is associated to the ver- 
tical mass flux: 

2 

6f = vJvz-\ -c6c, (8) 
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The perturbations of velocity, sound speed, density, pres- 
sure, vorticity are directly related to 5f, Sh, SS through a 
set of equations recalled in Appendix A (Eqs. (jAUIATp V 
The differential system is a simplified version of Eqs. (22- 
24) of FSJ06: 
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where the frequency Wev, related to the evanescent or 
propagating nature of acoustic waves in the vertical di- 
rection, is defined by: 

cjev = fcxc(l -M^)i (13) 

The boundary conditions at the shock surface are de- 
duced from Eqs. (28-30) of FSJ06: 
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where the velocity of the shock is related to its dis- 
placement through Av = —iujAC. The expression of 
the leaking lower boundary condition is deduced from 
Eq. (34) in FSJ06, which is based on the decomposition 
of a perturbation onto acoustic and advected contribu- 
tions (Eqs. (|M)) - (|XTT|) in Appendix A): 
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Fig. 2. — Eigenfrequencies of the most unstable modes for 
< rii, < 7 in a flow where /H = 0.1. The value of the cor- 
responding wavenumber is indicated next to each eigenfrequency. 
Positive values of u!i correspond to growth. 



The parameter /x is defined by a square root of a com- 
plex number. The choice of its sign is set by the require- 
ment that the amplitude of acoustic waves oc exp{ikfz) is 
bounded in the direction of their propagation (Eq. (C19) 
of FSJ06): 
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2.3. Illustrative examples as possible benchmark tests 

The boundary value problem is solved by integrat- 
ing numerically the differential system ()10m2p from the 
boundary conditions ()14m6p at the shock to the lower 
boundary, using a Newton-Raphson algorithm to find 
the discrete set of eigenfrequencies such that the lower 
boundary condition pT|) is satisfied. Fig.[2]illustrates the 
typical irregular shape of the eigenspectrum of this type 
of flows. The eigenfrequencies that are plotted corre- 
spond to the most unstable harmonic for each transverse 
number < ria; < 7. The following set of parameters has 
been used: cf^^/cl^t — 0-75, = 4, Hy = 0.1, where 
lengthscales are expressed in units of H. Eigenfrequen- 
cies are shown using a reference timescale Taac defined 
as the sum of the advective and longitudinal acoustic 
timescales between the shock and the deceleration re- 
gion: 
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Although the mode n^; = 1 is the most unstable in Fig.O 
this is not a strict rule. Fig. [3] shows that the order 
Ux of the most unstable mode depends on the size of 
the coupling region. For example, the linear instability 
should be dominated by a mode ria, = 4 if ffy = 0.04, 
and by Ux = 7 if i?v = 0.02. 
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Fig. 3. — Dependence of the growth rate uii on the size Hs/ of 
the coupling region, when the ampUtude A$ of the potential jump 
is kept constant. Same parameters as in Fig. [2] 

According to Fig. [3l low order modes seem to be 
favoured by the large size of the coupling region. This 
observation will be confirmed by physical arguments in 
Sect. 5. Note that the stability of the longitudinal mode 
rij, = is not an intrinsic feature of this toy model, but 
the consequence of tuning the size of the potential jump. 
In the limit H-^ —>■ 0, the mode 71^ = would become 
unstable if cf^/c^^t ^ 0.7167, but even then it would still 
be less unstable than transverse modes. Even for a weak 
potential jump such as cf^/cl^t = 0.9, some modes are 
still unstable (e.g. = 18) if the coupling region is 
sufficiently narrow. No unstable modes where found for 
cl/cl^^ > 0.95. 

3. METHOD FOR THE DETERMINATION OF THE 
INSTABILITY MECHANISM 

3.1. Cycles efficiencies Q and TZ 

The determination of the eigenfrequencies in Fig. [5] is 
not sufhcient to reveal the instability mechanism at work. 
Our toy model has been designed to be simple enough to 
allow for an unambiguous determination of the mecha- 
nism, using the same method as in Foglizzo (2002) for 
black hole accretion or Foglizzo et al. (2007) and YF08 
for neutron star accretion. In these studies, this method 
was restricted to perturbations with a high enough fre- 
quency for the validity of the WKB approximation in 
the vicinity of the shock. This condition excluded the 
description of waves propagating horizontally along the 
shock. By contrast, the simplicity of the present set up 
allows for an exact calculation even at low frequency. 
It leads to an analytical formulation of the advcctivc- 
acoustic instability in a compact approximation in the 
same spirit as Marble & Candel (1977) or Foglizzo & Tag- 
ger (2000) in ID. The originality of the present formula- 
tion is that it includes transverse motions and vorticity 
perturbations, since our model is 2D. This toy model 
can thus describe the vortical-acoustic cycle as well as 
the entropic-acoustic cycle. 



Our method consists in analyzing the interaction of plane 
waves (i.e. with a real frequency) with both the shock 
and the region of deceleration. Perturbations with a 
given frequency ujr and transverse wavenumber kx are lo- 
cally projected onto pressure waves and advected waves 
(entropy/vorticity) (see Appendix A). These different 
waves are independent from each other in each uniform 
part of the flow, but get coupled together both at the 
shock and in the region of flow gradients. This coupling 
is measured by two efficiencies T^sh and T^v of acoustic 
reflection, and two efficiencies Qsh, and Qv of advective- 
acoustic coupling, illustrated in Fig.[T] Connecting these 
waves into two global cycles, one can estimate the global 
efficiency Q = Qsh x Qv of the advective-acoustic cycle, 
and the global efficiency TZ = TZsh x TZ-v of the purely 
acoustic cycle in order to address the question of the sta- 
bility of each of these cycles. 

As illustrated by Fig. [U the coefficient Qy is the net ef- 
fect of the advective-acoustic coupling associated to the 
region of deceleration, together with the advection from 
the shock to this region and the acoustic propagation 
from it up to the shock. Similarly, the coefficient TZ^ 
includes the acoustic propagation between the shock and 
the region of deceleration, in addition to the acoustic re- 
flection by the deceleration region. 

Immediately after the shock, the perturbation of en- 
ergy density is decomposed into acoustic ones propagat- 
ing downward {5f^), upward {5f~^ and advected ones 

Sfsh = 6fi + 6f++Sf-^. (21) 

The definitions of the coupling constants TZsh, Qsh, TZ^, 
Qv correspond to the following relations: 

<5/+=7^.h<^/-„ (22) 
Sfi = Qsh6f^,, (23) 
Sf^,^QvSfl + TZ^Sf+. (24) 

Suppressing 5f^^, 5f^^ from these relations, 

Q + 7^ = l, (25) 

with 

Q^QshQv, (26) 
7^=7^,h7^v• (27) 

The analytic calculation of the coupling efficiencies Qsh, 
TZsh at the shock was done approximately in Foglizzo et 
al. (2005), and the exact calculation is repeated below 
in Sect. 4. 

The calculation of the coupling efficiencies Qv, TZv re- 
quires in principle an integration across the region of cou- 
pling, as was done in Foglizzo (2002). Taking advantage 
of the simplicity of the present flow, we show in Sect. 5 
that an algebraic expression of the coupling efficiencies 
Qv and TZ\/ is possible in the compact approximation. 
Anticipating on the results of Sect. 6, this calculation 
concludes that the purely acoustic cycle is always stable 
(i.e.|72.| < 1), and that only the advective-acoustic cy- 
cle may be unstable (i.e. |Q| > 1). This approach also 
proves that transverse perturbations are more unstable 
than longitudinal ones. 
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Fig. 4. — Efficiencies |Qahl a-nd [TZahl of the adve ctiv e-aco usti c 
coupling and acoustic reflection, deduced from Eqs. II29I I and I I31I I. 
Tliese efficiencies depend on the frequency of the perturbation only 
through the parameter /i^j^. 



3.2. Cycle of waves and growth rate: a new method to 
estimate the contribution of each cycle 

In Foglizzo (2002) and Foglizzo et al. (2007, here- 
after FGSJ07), the efficiencies \Q\, \TZ\ of the cycles of 
plane waves have been computed for a real frequency 
to — LUr- Eq. (|25|) invites us to evaluate the contribu- 
tion of each cycle considered alone through the complex 
frequency cu^'^ solution of Q{uj) — 1 for the advective- 
acoustic cycle, or the complex frequency u;^'^'^ solution 
of TZ{uj) = 1 for the purely acoustic cycle. The eigen- 
frequencies ui^'^'^^ 0;?^"= associated in this manner to each 
cycle can be compared to the eigenfrequency lo computed 
directly in Sect. 2, in order to assess their relative impor- 
tance in the instability mechanism. This method is used 
for the first time in Sect. 16.2] 



4. EFFICIENCY OF ACOUSTIC REFLECTION AND 
ADVECTIVE-ACOUSTIC COUPLING AT THE SHOCK 

4.1. Acoustic reflection at a shock 

The efficiency of acoustic reflection has been computed 



by Foglizzo et al. (2005) (Eq. (FIO)): 
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The WKB approximation in Eq. (|B5|1 of Appendix B 
(or Eq. (CI) in Foglizzo 2002) indicates that p,\Sf^\'^/vc 
is conserved to second order in an adiabatic flow. This 
guarantees that the quantity |7?.shP can also be inter- 
preted as a ratio of fluxes of acoustic energies. 
From the expression of kf in Eq. (jl9p , acoustic waves are 
propagating in the z direction if /i^ > (i.e. to > lUcv) 
and evanescent otherwise. Using Eq. ((29|) . one can prove 
that \TZsh\ < 1 for propagating acoustic waves (^^j^ > 0) 
and |7?.sh| = 1 for evanescent ones (^^j^ < 0), as illus- 



trated by the dotted line in Fig. |4l This property pre- 
cludes the possibility of an over-reflection, at least in this 
simple toy model. 

4.2. Advective-acoustic coupling at a shock 

The coupling coefficient obtained in the WKB approx- 
imation in Foglizzo et al. (2005) (Eq. (F11-F12)) can be 
rewritten as an exact formulae in the present toy model 
since the flow is uniform immediately after the shock: 
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This formula is checked by numerical simulations in pa- 
per II. Note that |QshP is not a ratio of energy fluxes 
(unlike |7^shP). 

Eq. PT|) indicates that the efficiency of the advective- 
acoustic coupling vanishes in the limit of a weak shock: 
Qsh oc (1 — Algij) for radial perturbations (/ish — 1), and 
Qsh oc (1 — A^sh)^ fo"^ transverse ones (/Ush < !)■ As seen 
on Fig. m there are two optimal values of cx ±1/A(^ 
leading to a maximal production | Qsh | of advected per- 
turbations. In the asymptotic limit of a strong shock, 
the maximum efficiencies are reached for acoustic waves 
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5. COMPUTATION OF THE ACOUSTIC FEEDBACK 

5.1. Analytic expressions ofTZ^, Qv *^ the compact 
approximation 

5.1.1. Range of validity 

The "compact approximation" consists in treating the 
inhomogeneous layer as infinitely thin. This approxima- 
tion was used by Marble & Candle (1977) in the context 
of jet nozzles, in order to compute the acoustic feed- 
back from advected entropy perturbations, in ID. It is 
here extended to parallel flows in an external potential 
$, subject to both entropy and vorticity perturbations. 
This approximation is valid for perturbations such that 
their vertical wavelength 2n/kz is long compared to the 
lengthscale Hy of the region of deceleration, on both 
sides of it. For acoustic waves, the criterion k^H^ <^ 2n 
can be approximated, using Eq. (|19|) and Cm < Cout, by 
the following sufficient condition on the frequency lu of 
the wave: 
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For advected waves {k^ = oj/v), the frequency threshold 
for the compact approximation can be related to the time 
Tv of advection through the inhomogeneous region: 



27r 



> TV 



dz 



(34) 



6 



1.2 



0.8 - 



0.6 



0.4 



0.2 



1 . . . 

IR 1 (H /H=0) 

V V 

- IR 1 (H /H=0.1) 

V V 

- n =1 

X 1 

- n =0 

X 




n =5 

X 


: i 

f 


■ # 

■ f 
• f 

i 


/\ ■ 



H 
0.01 



1.2 



0.1 



10 



0) 



Fig. 5. — Refraction efficiency |7^vl foi' longitudinal and trans- 
verse perturbations, in a flow where H^/H = 0.1. The differ- 
ent curves are labelled with the value of the transverse wavenum- 
ber Tlx = 0,1,5. The thick grey lines are the analy tica l esti- 
mates deduced from the compact approximation, in Eq. (|35^ . The 
range of frequen cies corresponding to total reflection coincides with 
[4^,a;->'] (Eq. GSj) 

The condition p4p for advected waves is more severe 
than ([55)1 because the flow is subsonic. The coefiicients of 
acoustic refraction, and advective-acoustic couphng can 
be expressed through simple algebraic formulae based on 
the conservation of entropy, mass flux and energy density 
{Sf, Sh and SS) through the potential jump. 

5.1.2. Coupling ejficiencies TZv and Qv 

In Appendix C, the following expressions are deduced 
from the conservation laws : 
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where tq and t-ji are defined by: 
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As expected, these coupling efficiencies vanish in the 
limit of a uniform flow (cout — Qn and Eq. ([5])). 
The exponential functions in Eqs. (|35ll36p are associated 
to the vertical structure of advected and acoustic waves. 
In the low frequency limit (w <C w™), the evanescent 
character of acoustic waves is responsible for the follow- 
ing damping effect: 
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Fig. 6. — Efflciency |Qvl of the advective-acoustic coupling in 
the same flow as in Fig. [5] The thick grey lines corresponds to 
Eq. II36I I. Kinks in the curves correspond to the particular frequen- 
cies lijj," and . There is no frequency cut-off in the compact 
approximation. 
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This acoustic evanescence is a strong stabilizing factor 
for any cycle based on high degree perturbations {rix ^ 
L^/2TrH). 

5.2. Comparison with the exact calculation: the 
frequency cut-off uj^ 

Figures [5] and [6] show the efficiencies \Qv\, \T^v\ com- 
puted in the compact approximation (thick grey lines) 
and the exact calculation with H\j/H = 0.1 (thin black 
lines). The numerical method to compute Qvi is 
based on three numerical integration from the shock to 
the lower boundary, using different boundary conditions 
at the shock, and two appropriate linear combinations 
of these three solutions such that the lower boundary 
condition is fulfilled, as explained in Appendix D of 
FGSJ07. Qv can also be computed as an explicit integral 
over the region of flow gradients, using Green functions 
as explained in the Appendix D of the present paper. 
Schematically, the global advective-acoustic coupling Qv 
is described by a convolution between the radial profile of 
acoustic waves 5pQ /p, the radial profile of advected waves 
q! ^dz^ local emissivity dh\j /dz associated to the 

gradients of the flow: 



Qv= / oo e^^h " -7:—a-z 

Jhc P oz 



(41) 



where 
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, lU) 
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dz 
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(43) 



In comparison with Eq. (25) of Foglizzo (2001), the 
present formulation stresses the focal character of the 
acoustic emissivity associated to the flow gradients, de- 
scribed by Eq. (gS]). 

Figures O and [B] show that the analytical compact ap- 
proximation is excellent for cjtv <C 1, and works quite 
well even for values just below one. As expected, r^^ 
defines a cut-off frequency above which the compact ap- 
proximation ceases to be valid. From Fig. [6l one can 
estimate the cut-off associated to the advective-acoustic 
coupling: 



1 

TV ' 



(44) 



Comparing Figs. [5] to Fig. [HI the frequency cut-off as- 
sociated to Q is smaller than the one associated to 7?., 
as expected from Eqs. (pSj) and pi]) . The different na- 
ture of these cut-off frequencies can be understood owing 
to Eq. (j4T|) : the efficiency Qv decreases not only if the 
wavelength of the perturbations is short compared to the 
lengthscale of the flow gradients (i.e. the term db^ /dz), 
but also if the convolution of the advected and acous- 
tic phase structures leads to a phase averaging. Also 
remarkable in Fig. [5] is the fact that acoustic waves prop- 
agating horizontally {uj ~ wj") are always well described 
by the compact approximation, with |7^v| = 1 (this is 
also visible in the bottom plot of Fig. [7]). Altogether, 
Figs. [5] and [H] indicate that the compact approximation 
may be used as an upper bound for the efficiency of the 
acoustic feedback. 

Acoustic reflection is always damped (|7^v| < !)• Total 
reflection without damping (|7^v| — 1) occurs for trans- 
verse perturbations in the range of frequencies [wj," , 
where both /i?^ > and /io^j < 0, i.e. when acoustic 
waves propagate above the inhomogeneous region and 
are evanescent below it. 

A striking feature on Fig. [B]is the divergence | Qy | oc 
of the coupling efficiency at low frequency for transverse 
perturbations (Eq. ([36])). This divergence however is 
not sufficient to overcome the weak efficiency of produc- 
tion of advected perturbations at the shock, |Qsh| oc 
(Eq. (pij) ). The frequency dependence of the global ef- 
ficiency |Q| = IQshQvl is further studied in the next 
section (Fig. [T]). 

6. GLOBAL CYCLE 

6.1. Compact approximation of Q and TZ 

The coefficients Q and TZ defined in Eqs. (|26ll27p can be 
expressed analytically by combining Eqs. ([31]) and (|36p 
for Q, and Eqs. (EH) and ^ for TZ: 



-TZoe''^^'^, 



(45) 
(46) 



Mil 



1 



7+1 Min (l-Mi„Xi„)(M?n 
A^out+Mout {1~M?J{1^ 



Mout 
1 + Mout A^out 



f 2/iinXi„-)-Ali 



Mout" 



Ma 



(47) 
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Fig. 7. — EfScicncics \Q\ and \Ti\ of the advective-acoustic cycle 
in a flow where /H = 0.1 (black lines). The different curves are 
labelled with the value of the transverse wavenumber < n^: < 7. 
In the upper plot, the crosses indicate the discret set of oscillation 
frequencies corresponding to the solution of Q = 1 (Fig. [9]l. The 
grey lines are the analytical estimates deduced from the compact 
approximation. For comparison with Fig. (6] Taac/Tv ~ 16.52. 



TZo = - 



1 - flinMin) fif^ + 2^i„Xin + 



Mout A^inC? 



MinA^outCout + MoutXinCf„ 



(48) 



The cycles efficiencies |Q| and \TZ\ are shown as functions 
of the real frequency tOr in Fig. [7l both in the compact 
approximation and for H^/H = 0.1, for different values 
of the transverse wavenumber n^- The stability of the 
purely acoustic cycle is confirmed by the bottom plot in 
Fig. [71 The spikes correspond to the frequencies wj,", for 
which \TZ\ = 1 according to Eqs. ^ and ([iSjl . 

The stabilization of the advective-acoustic cycle above 
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Fig. 8. — Eigenfrequencies of the most unstable modes for 
< n^; < 7 in a flow where H\; = (empty diamonds). The 
value rix of the corresponding wavenumber is indicated next to 
each eigenfrequency. Increasing the size i^v of the deceleration 
region up to jtl = 0.1 is stabilizing (filled diamonds). 
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Fig. 9. — For each value of = to 7, comparison of the most 
unstable modes associated to the advective-acoustic cycle alone 
(filled grey squares), the purely acoustic cycle alone (black crosses) 
and the full problem (filled black diamonds). The contribution of 
the purely acoustic cycle can either be constructive or destructive. 



the cut-off frequency is clearly demonstrated by the up- 
per plot in Fig. [T] The maximum coupling efficiency is 
always reached for w w^", as expected from Sects. 4 
and 5 since usf^ — uf^ in our toy model. If the shock is 
strong and A^out 1, we deduce from Eq. ([17)) that the 
maximum efficiency is 

Qmax^^ -i (49) 

7+1 XL 4 

In our toy model with a strong shock, the maximum 
efficiency Qmax is larger in 2D than in ID by a factor 
which depends only on M.\i^: 

gmax(2J) ^ 1 + 2.Mi„ 

This factor is equal to 1.8 for 7 — 4/3. 
According to Fig. [7l the advective-acoustic cycle alone is 
necessarily stable (|Qmax| < 1) for n^; > 4 if H\j / H = 
0.1, and can be unstable for rix < 3. Interestingly, Fig. [5] 
indicates that the mode n^; = 2 is stable and the mode 
rij; = 4 is unstable, whereas the upper plot of Fig.[7]shows 
that \Q\ > 1 is possible for the mode Ux — 2, and that 
|Q| < 1 for Ux = 4:. This apparent paradox is clarified 
below by taking into account two factors: 

(i) the phase closure relation (|25p selects a discrete 
set of oscillation frequencies, which do not necessarily 
allow for the maximum coupling efficiency |Q|max to be 
reached, 

(ii) the additional influence of the purely acoustic 
cycle can be constructive or destructive, and can be 
crucial close to marginal stability. 

The closure relation ([^5]) takes the following particular 
form in the compact approximation: 

Qoe*""« -f Teoe*""^ = 1, (51) 



This explicit analytical equation is a 2D generalization of 
the equation obtained in ID by Foglizzo & Tagger (2000). 
The corresponding eigenmodes, easily obtained numeri- 
cally, are shown as empty diamonds in Fig. [51 Comparing 
the compact eigenspectrum (filled diamonds) to the one 
obtained for H-^ / H = 0.1 (empty diamonds) shows again 
the stabilization above the cut-off frequency. 

As a consistency check, one can verify for 1 < Ux < 1 
that the growth rates shown in Fig. [31 in the limit 
H^/H ^ 1, do converge towards the growth rates rep- 
resented as empty diamonds in Fig. [51 The former are 
obtained by numerical integration of the boundary value 
problem, whereas the latter are directly deduced from 
the exphcit Eqs. (|47l48l5ip . 

6.2. Calculation oj uj'^'"' , uj"^"^ 

As explained in Sect. 13.21 the complex eigenspectrum 
associated to each cycle considered alone is defined by 
the equations Q = 1 for the advective-acoustic cycle, 
and TL = \ for the purely acoustic cycle. These eigen- 
spectra are compared to the solution of the full prob- 
lem {Q + Ti = 1) in Fig. [31 where we have selected 
the most unstable mode for each degree < n^: < 7. 
The eigenspectrum of the purely acoustic cycle is stable 
as expected. The eigenmodes of the advective-acoustic 
cycle match remarkably well the eigenmodes of the full 
problem from the point of view of the oscillation fre- 
quencies, and also the growth rates for < n^; < 4. 
The difference of growth rate between the solutions of 
Q + TZ — 1 and Q = 1 is attributed to the contribution 
of the purely acoustic cycle, which is clearly constructive 
for the modes Ux — 0,1,4,6,7 and destructive for the 
modes Ux — 2,5. As already noted by Foglizzo & Tag- 
ger (2000) and Foghzzo (2002), the contribution of the 
purely acoustic cycle is most important close to marginal 
stability: Fig. [51 shows that the mode Ux = 2 would be 
slightly unstable with the advective-acoustic cycle alone. 
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Fig. 10.— Same as Fig.[3l with L^/H = 3.5 instead of L^/H = 4. 
This small change in the aspect ratio of the flow changed the degree 
Ux of the dominant mode: the modes Jia; = 1, 2 and 5 are now the 
dominant ones depending on the size of the region of coupling. 



and is stable because of the destructive interference with 
the acoustic cycle. Conversely, the mode n^^ = 4 would 
be stable with the advective-acoustic cycle alone, and 
is unstable because of the constructive interference with 
the acoustic cycle. 

Unlike the method based on Q{LOr) (e.g. Foglizzo 2002, 
FGSJ07), the present method is not restricted to uoi <C 

6.3. Analytical estimate of the maximum growth rate 
wf'^'^ and sensitivity to geometrical factors 

Writing oj'^'^'' = uj^^'^'^+iui'^'^'^ , the growth rate wf'^'' of the 
advective-acoustic cycle alone can be estimated by a first 
order expansion of Qo exp(a;'*'^'^TQ) = 1 in the compact 
approximation: 



^lr + Mi, 



— log|Qo« 



Mr 



;1 



(52) 



(53) 



This equation is valid when ujf'^'^ <^ w^'"^ and > 0. 
Equation (|52p illustrates the competing effects intro- 
duced by transverse perturbations. On the one hand, 2D 
perturbations allow for a higher coupling efficiency Qq 
(Fig. [7] and Eq. [5D]), but on the other hand the acoustic 
timescale for the feedback is longer for transverse waves 
(factor Hr in Eq. ([5^ ). Whether the optimum balance 
between these two effects is reached or not depends on 
geometrical factors such as the aspect ratio Lx/H of the 
toy model. Indeed, the discrete set of oscillation fre- 
quencies selected by the closure relation is sensitive to 
Lx through kx = 2'Knx/Lx- Besides, the constructive or 
destructive contribution of the acoustic cycle is sensitive 
to H through the ratio of the acoustic and advective- 
acoustic timcscales. A comparison between Fig. [3] and 
Fig. [To] illustrates the fact that a small change in the 



aspect ratio of the unstable cavity, from Lx/H = 4 to 
Lxl H — 3.5, is sufficient to favour the modes Ux — 1,2,5 
rather than Ux = 1, 4, 7. 

6.4. A low frequency, low order instability if the 
coupling region is wide 

The frequency cut-off associated in Sect. l5.2l to the 
finite size of the coupling region is responsible for the sta- 
bilization of high frequency modes. Besides, the evanes- 
cent pressure feedback associated to high Ux transverse 
modes (Eq. |40)) is inefhcient at perturbing the shock. 
These two arguments explain why the advective-acoustic 
instability favours low frequency, low order modes if the 
size of the coupling region is large, as seen on Figs. [Hand 

m 

More precisely, the evanescence length \ha{k^ )] ^ of 
the pressure perturbation can be compared to the dis- 
tance H to the shock. Using Eq. p9)) . the conditions 
Im(fc^i/) < 1 and lj < t^y can be translated into a con- 
dition on the horizontal wavenumber kx : 



ki < 



1 - 



(54) 



Neglecting A^?^ <C 1 for the sake of simplicity, and us- 
ing Eq. (|44p . the range of frequencies and wavenumbers 
which can produce an unstable cycle most efhciently are 
summarized as follows: 



1 



nl < 



TV 

Lx 



27rCinTV 



+ 



/ Lx 



\2'kH 



(55) 



(56) 



The r.h.s in Eq. (|55)) depends both on the velocity profile 
in the deceleration region (phase averaging), and on the 
shock distance (acoustic evanescence). The first term is 
dominant in the compact approximation. 

7. RELATION TO THE INSTABILITY OF THE STANDING 
SHOCK DURING CORE COLLAPSE 

A direct extrapolation of Eq. (|56p to a stalled shock of 
radius i?sh in spherical geometry would favour unstable 
modes with a low azimuthal degree /: 



Rsh 

H 



(57) 



The deceleration region responsible for the acoustic feed- 
back has been identified in the numerical simulations of 
Scheck et al. (2008), from which we estimate, 80ms after 
bounce in their model WOOF: Rgh ^ 110km, H ^ 50km, 
-Rsh/cin ~ 5ms, Ty 15ms and / H ^ 0.4. Using 
these numbers in Eq. (|57)l we deduce I < 2.2. Despite 
the simplicity of our approximation, this estimate is very 
close to the degree I — 1,2 of SASI observed in Scheck 
et al. (2008), as well as in most numerical simulations of 
the stationary shock involved in the core-collapse prob- 
lem. These simulations cover a broad range of physical 
ingredients and approximations, from stationary flows 
with a simplified cooling function (Blondin & Mezza- 
cappa 2006, FGSJ07) to the most detailed modeling of 
neutrino transport and advanced equations of state in the 
core of a massive progenitor (e.g. Burrows et al. 2007, 
Marek & Janka 2007). SASI has also been observed in 
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the pioneering 3D simulations of Blondin & Mezzacappa 
(2007) and Iwakami et al. (2008). The efficiencies Q, 
TZ of the cycles have been computed by FGSJ07 in the 
context of a simplified stationary flow, proving that the 
advective-acoustic mechanism is responsible for SASI, at 
least when the shock is far enough to allow for the WKB 
approximation. By taking a step further in the simplifi- 
cation of the flow, the present toy model enabled us to 
better understand some generic properties of this insta- 
bility. A frequency cut-off is associated to the timescale 
of advection through the flow gradients. Among the 
low frequency perturbations below this cut-off, the most 
unstable ones are those which approach the maximum 
advective-acoustic coupling efficiency Q with a moder- 
ate lengthening of the acoustic feedback timescale. Close 
to marginal stability, constructive interference with the 
purely acoustic cycle can be essential for the instability. 
In such cases however, its consequences may be negligi- 
ble due to a lack of time for significant growth. 
By considering a 2D planar flow, our toy model neglected 
the large scale gradients that would exist in a spherical 
flow, all the way between the shock and the neutron star. 
These large scale gradients are expected to contribute 
the adiabatic coupling of advected and acoustic pertur- 
bations at low frequency. Neglecting these gradients in 
the toy model was important in order to characterize 
analytically acoustic waves in the vicinity of the shock, 
even at low frequency (Eqs. (jAlOHAlTj) ). Another impor- 
tant simplification in our toy model was the hypothesis 
of an adiabatic region of deceleration, which allowed us 
to identify the phase averaging effect responsible for the 
inefficiency of the acoustic feedback, above the frequency 
cut-off, through Eq. (|4T|) . 

Of course, the acoustic feedback produced by advected 
perturbations is expected to depend quantitatively on 
whether the deceleration is produced by a cooling layer 
or an adiabatic gravity step. For example, the mode 
I = would be necessarily stable in an idealized isother- 
mal flow, because the only advected perturbations which 
can build a cycle in an isothermal flow are non radial vor- 
ticity perturbations {I > 1). The relative contributions of 
entropy and vorticity perturbations to the acoustic feed- 
back are thus expected to differ significantly from our 
toy model when non adiabatic processes are involved. 
We believe however that our description of the cycles, 
and the existence of a cut-off frequency associated to the 
advection time through the deceleration layer, are gen- 
eral enough to be relevant to the framework of stalled 
accretion shock during core collapse. This toy model al- 
lows us to propose answers to the questions: why SASI 
is a low frequency, low degree I instability, and why the 
eigenspectrum is irregular and sensitive to geometrical 
factors. 

8. CONCLUSION 



A simple toy model has been analyzed in order to bet- 
ter understand the advective-acoustic cycle in a decel- 
erated flow. This comprehensive study has reached the 
following conclusions: 

-The instability mechanism is based on the amplifica- 
tion of the advective-acoustic cycle. 

-The purely acoustic cycle is always stable in this flow. 

-Optimal transverse perturbations are more unstable 
than longitudinal ones. 

-The finite size of the coupling region is responsible for 
a frequency cut-off which favours low frequency pertur- 
bations. 

-This frequency cut-off disfavours high order perturba- 
tions, which would lead to an inefficient evanescent pres- 
sure feedback: the most unstable modes are expected to 
be low frequency and low order. 

-The maximum efficiency |Q| of the advective-acoustic 
cycle is governed by a compromise between the efficiency 
of the advective-acoustic coupling, which favours trans- 
verse perturbations, and a short cycle time which dis- 
favours too horizontal acoustic propagation. 

-The eigenspectrum is irregular and sensitive to geo- 
metrical factors due to the phase condition which selects 
a discrete set of frequencies, and the interferences with 
the acoustic cycle. 

-A new method has been described to measure the 
eigenspectrum associated to each cycle considered alone. 

-This toy model can be used as a benchmark test for 
numerical simulations, in order to check for undesired 
numerical artifacts. The important roles played by the 
shock, the advection of perturbations and the propaga- 
tion of acoustic waves make this toy model particularly 
challenging for numerical simulations (paper II). 

-The comparison with the properties of SASI suggests 
that our toy model captures some fundamental features, 
explaining in particular why SASI is a low frequency, low 
degree instability. 

A good understanding of the mechanism behind SASI 
should help perform accurate numerical simulations of 
core collapse explosions, with particular attention to the 
mesh size and the boundary conditions. 
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APPENDIX 

DECOMPOSITION OF PERTURBATIONS IN A UNIFORM FLOW 

For completeness we repeat here the relationship between Sf, 5h and Sv^, Sp, 5p corresponding to Eqs. (A1-A4) in 
FGSJ07: 
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Sp 



IP 



(A3) 



^ = (A5) 

The expression of the perturbations of transverse velocity and vorticity are simply proportional to Sf and SS, as 
deduced from Eqs. (A15-A18) of YF08: 

Sv^^^Sf, (A6) 

Swy = '^ -5S. (A7) 

w 7 

In a uniform flow, the perturbations of energy density 5f^ and mass flux 5h^ associated to advected perturbations of 
entropy /vorticity produced by the shock are given by Eqs. (C5-C6) in FGSJ07: 

c2 1 - fi-^M^ 7 ' ^ ' 

Sh'^ = i^Sf^ - SS. (A9) 
The energy density Sf^ and mass flux Sh^ in the acoustic waves are deduced from Eqs. (C7-C8) of FGSJ07: 

6f^ = b±^{6h + 6S)- 6f, ( AlO) 



2 2/i 

Sh^ = ±i.'-^. (All 

The perturbations of pressure Sp^ and energy density Sf^ of an acoustic wave in a uniform flow are related by: 

Sp± Sf^ lT^^M 



7P 1 — 

WKB APPROXIMATION 



(A12) 



Owing to the adiabatic hypothesis, the vertical structures of both entropy and vorticity perturbations can be inte- 
grated and the differential system can be reduced to a second order differential equation, similar to Foglizzo (2001) 
when considering the perturbation of energy density Sf and the variable X defined by: 



dX 



dz l-7V/(2' 



(Bl) 



„ 1 dSf 5f f 1 \ 5S 

iLu oX & \Ai / 7 



This equation is formally identical to Eq. (B32) of Foglizzo (2001) (with 6K ~ 0). This formulation is particularly 
useful to define a WKB approximation Sf^y.^^ of the acoustic solutions of the homogeneous equation in the simplest 
manner, satisfying a WKB criterion deduced from the left hand side of Eq. (jBSP : 



^»il-M^f-^(4-,], (B4) 
c oz yv^c^ ' 

Sh^^±^^J±, (B6) 

where Sf^ is an arbitrary normalization constant. The acoustic solution satisfying the lower boundary condition 
defines the reflection coefficient TZ-^, such that Sfo = 6f + + TZ^Sf^ at z = z^h- 
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COMPACT APPROXIMATION 
Dispersion relation 

The dispersion relation is obtained as in Foglizzo & Tagger (2000), Foglizzo (2002). In the specific case of our 
toy model in the compact approximation, some coupling constants , TZ'-^ can be defined immediately above the 
coupling layer (2; ~ zy) by 

Sf^ ^QvSM + n'^Sf+. (ci) 

The value of perturbations at 2; = Zsh and z = zy are related through their vertical wavenumber kf defined by Eq. (|19p 

and fc^ = Lu/v. 

Sfsi.^Sfi + Sf++5f-, (C2) 
Sfi = e^'""Sfl (C3) 

^f+^^rktH^j+^ (C4) 

6f-^^e^KHsf^_ (C5) 

The coupling constants Qy and T^y can thus be decomposed into a contribution of advection/propagation, and a 
contribution of the coupling through the compact layer: 

Qv-Q'ye*""«, (C6) 

7^y=7^Ve''^^^, (C7) 

with 

TQ = [k, ~ fcj)— = Taac , , (C8) 

Note that the quantities tq and t-h may be directly interpreted as timescales only when the frequency uj is real and 
high enough to correspond to propagating acoustic waves {fif^ > 0). The global efficiencies Qo and TZq are defined by 
Qo = QshQ'y and 7^o = 7^sh7^V■ 



Reflection of acoustic waves: calculation of T^y 

Let an acoustic wave propagate in the direction of the flow towards the region of gradient. It produces an acoustic 
reflection Sf^^ and a transmitted acoustic wave 

6fiL + Sfr = Sf+„ (CIO) 
Sh++6hr^^5ht^,. (Cll) 

The system (jClOtiClTj) is thus equivalent to 

Sf^. + Sfr,^Sf+,., (C12) 



The reflexion coefficient TZy in Eq. ([55]) is deduced. 



^/ii - sf,- = T ^fouf (C13) 



Advective- acoustic coupling: calculation of Qy 

At low frequency, the acoustic perturbations produced by the advection of the entropy /vorticity perturbation through 
the gradient satisfy the conservation of mass and energy: 

Sffn+Sfrn=S.fL+SfL, (C14) 

5hi + dK^^Sh^^, + Sht^,. (C15) 



Using Eq. (Km . 



Sfin-Sf^^.^Sfi.^Sfg, (C16) 

/^in~I^ ~ Mout 2°"* • (C^l''') 



''out 



thus, the advective- acoustic coupfing Qy is described by Eq. ([55)1 . (5/o^t is deduced from 6f^^ by exchanging the 
subscripts 'in' and 'out', /iin/A^in and — /iout/A^out- 
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LOCAL ACOUSTIC EMISSIVITY 



Expressing the local pressure emissivity 05 characterizing to the advectic-acoustic coupling due to the gradients 
of the flow is most easily obtained by writing the differential equation satisfied by pressure perturbations Sp. A 
pressure perturbation in an adiabatic flow satisfies the acoustic wave equation with a source term associated to the 
inhomogeneity of the flow. The differential equation established by Foglizzo (2001) in a spherical accelerated flow is 
formally the same in a decelerated flow. The differential system can be written as follows: 



ai 



d_ 

dz 



a2 



Sp 
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asSS, 



with 



A; 
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dlogM 
dz 
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dz 



iiov^ 1 + A^^ 



1 -7^2 



(Dl) 

(D2) 
(D3) 

(D4) 
(D5) 
(D6) 



Let be the acoustic solutions of the homogeneous equation which propagate outward (-) or inward (+), and let 
(5po = Spj^ +TZ-^dp^ be the acoustic solution which satisfles the lower boundary condition at rbc- The general solution 
satisfying the lower boundary condition is written as 

f as f, ^^, 1-M^ Spo^ ( 



5p 



P 



as f -dz 1 - Sp^ 



where W contains the normalization of the amplitudes of Sp± : 

\6po d fl-M^ 5p- 



p dz 



P 



Sp^ d 
p dz 



P 



X2 



5po 
P 



d}ogWp 
dz 



= -ai. 



The WKB approximation of the pressure perturbation deduced from Eq. (IBSP is 



Sp: 



wkb 



IP 



1-X2 c2 ' 



2 A'sh 



1 -A42 



Msh ujc. 



sh 



Thus 



as 



The derivative of the pressure perturbation is: 



izd_ 
dz 



c 
A 



2v 



d Sp _ d ^ 
dz SSsh dz 



as r, . 

pJsh 



'-dz 



7W2 



—dz] - -^Spo [ f3 
p J dz ' 



dlogM 
dz 



'-dz 



Sp^ 

A^2 p 



The outgoing pressure perturbation is deduced from the WKB approximation 

Sp — Sp^ + Sp^ , 

M 



d - iu 
—dp= — 
dz c 



Sp+ 



Each component Sp is deduced from Sp and its derivative according to: 



Sp± = 



UU±^)sp^-!-il^M-)^Sp\ 
2 y\ fi J lujfi dz J 



(D7) 

(D8) 
(D9) 

(DIO) 
(Dll) 
(D12) 

(D13) 

(D14) 

(D15) 
(D16) 

(D17) 
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The value of f3 is deduced from the outer boundary condition Sp- = 0, and the general solution satisfying the boundary 
conditions is thus simply: 



Sp 

ssT 



sh 



~ Sp- 



as f a^dz Spo 
p 



as r, ^dz\-J^ ^P- 

pJsh V 



■dz. 



The acoustic feedback produced by the advection of the entropy-vorticity perturbation is 

Spsh 



Qv = (1 - A*shA^sh) 



pSSsh ' 

1 -7^v 



■x ''sh 



p 



1 + 7^v 



from which we deduce Eqs. (|4HI43p . 



2d 9^ 



— w 



(D18) 



(D19) 



(D20) 
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